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FOREWORD 


This report presents the first year results of. an investigation being 
conducted under Contract NAS8-26363 for NASA George C. Marshall 
Space Flight Center under the technical direction of the Aero-Astro- 
dynamics Laboratory, Dynamics and Control Division. Dr; S. Winder 
is the technical monitor. The study is being performed by Convair 
Aerospace Division of General Dynamics under the direction of Mr. 

R. Huntington, project leader. 

The author is indebted to Messrs. R. Peloubet and R. Haller for their 
assistance in developing the computer program and in performing the 
numerical analysis and to Mr. J. Kramer for converting the program 
to run on the Univac 1108 computer. 
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SUMMARY 

A computer projf'ay*v~ wis been developed and demonstrated that can ana- 
lyze the response of space shuttle to atmospheric turbulence. The method 
developed accounts for propellant slosh, gimballed engine, and stability 
augmentation system (SAS) coupling with the elastic vehicle.. Statistical 
outputs are generated that relate vehicle loads and accelerations to level 
of random turbulence. For discrete turbulence descriptions, response 
time-histories are, computed. 

The response problem is formulated in the frequency domain, thereby 
permitting the use of any quasi-steady or unsteady aerodynamic theory 
based on harmonic motion. Gust penetration effects are readily accounted 
for. The quasi-steady approach using the Woodward aerodynamic method 
for subsonic and supersonic flows appears promising for the space shuttle 
turbulence response problem. Properly predicting flow interference, 
thickness, angle-of-attack, and body effects is more important for space 
shuttle than including aerodynamic lag effects for the reduced-frequency 
range of interest, ty 

Preliminary turbulence response analyses of space shuttle were conduc- 
ted for one ascent and one booster flyback subsonic flight condition. The 
ascent case considered both symmetric and antisymmetric boundary con- 
ditions, while only the symmetric analysis was conducted for booster fly- 
back. The three conditions were analyzed with SAS active and inactive, 
but apparent instabilities with SAS active (arising from improper gains) 
invalidated the results. Load and acceleration responses in all cases, 
however, were well within vehicle design limits. 

It is recommended that future response studies be preceeded by stability 
checks to ensure that the elastic vehicle/SAS is stable and that the gains 
used in the response analysis are correct. It is also recommended that 
the study be expanded to include other subsonic and supersonic flight 
conditions. 
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SECTION 1 


INTRODUCTION 

This report presents the results of the first year's work under Contract NAS8-26363, 
"Aeroelastic Effects on Shuttle Vehicle Dynamics". The objective of this portion of 
the effort was to develop and demonstrate a computer program to analyze the response 
of space shuttle to atmospheric turbulence. 

Section 2 contains the equations of motion, a discussion of unsteady and quasi-steady 
aerodynamic theories, and a description of methods for handling propellant slosh and 
. engine dynamics. The computer program is described briefly in Section 3; a detailed 
user's guide is contained in Reference 1. 

To check out and demonstrate the computer program, a typical space shuttle configura- 
tion was analyzed in the ascent and booster return conditions. The results are presented 
in Section 4. 

Application of random turbulence response techniques to launch vehicles has been very 
limited prior to the advent of space shuttle. For expendable systems in which fatigue 
life is not an important consideration and exposure to atmospheric turbulence is of 
sho,rt duration, discrete gust concepts have proven adequate for design. The space 
shuttle, however, will be subjected to many reuses and will encounter turbulence during 
ascent and flyback that may produce design loads, fatigue damage, and ride discomfort. 

Gust criteria for structural design have evolved from simple descriptions of atmospheric 
turbulence to more sophisticated concepts that correlate closely with measurements. 
Early design concepts treated the vertical (lateral) component of the gust velocity as a 
step function as shown in Figure 1-1. The aircraft was assumed to travel instantan- 
eously from a region with zero gust velocity into a constant gust velocity region. Sub- 
sequent refinements considered the effect of gradual penetration of the vehicle into the 
gust field. 

Later concepts considered one-dimensional discrete profiles such as the (1-cosine) 
wave shape shown in Figure 1-1. Efforts were then largely directed toward defining 
the wave length and peak velocity of the gust for aircraft design purposes. It was soon 
recognized that the average peak gust velocity associated with large wave lengths was 
greater than that for small wave lengths. 

The next step in the evolution of gust criteria treated the vertical (lateral) component 
of the gust velocity as being random in the flight direction while continuing to retain 
the concept of constant gust velocity in the direction perpendicular to the flight path 
(Figure 1-1). With the introduction of the random concept, gust velocity was defined 

. f 1-1 •’ 



by certain averages, such as root mean 
square (rms) and power spectral density 
(PSD). The observed fact that even the 
average values of turbulence vary with 
such parameters as atmospheric condition, 
geographic location, season of the year, 
and distance above the earth's surface was 
described in terms of the probability of 
encountering turbulence with specified rms 
values. By describing the aircraft re- 
sponse characteristics in terms of trans- 
fer functions, its response to one-dimen- 
sional random turbulence could be computed. 
The random response was also described 
in terms of averages, and it was possible 
to compute the average or expected number 
of times that the response exceeded any 
specified level. These concepts have been 
developed to the point that they are now 
widely accepted as the preferred manner 
in which to determine the effects of at- 
mospheric turbulence on aircraft. 

Numerous Air Force sponsored gust- 
measuring programs have been conducted 
and some are currently being performed. 
The results of these investigations have 
led to the generally accepted conclusion 
that the measured gust spectra can be 
best described by the Von Karman gust 
spectrum (compared in Figure 1-2 with 
the Dreyden spectrum) with a scale of turbulence (L) that varies from 500 feet near 
sea level to 2500 feet at higher altitudes. The probability of encountering rms gust 
velocities of any given value in clear air or thunderstorm turbulence is specified by 
the altitude-dependent variance of the rms gust velocities. The percentage of time 
that clear— air turbulence or thunderstorm turbulence is encountered is described by 
fractions, which also vary with altitude. 

The latest development of these concepts considers each component of turbulence to 
vary in any plane through the flight path of the vehicle. This concept, identified as 
two-dimensional random turbulence in Figure 1-1, is refinement of the one-dimensional 
turbulence model. It has been found, however, that the improvement in the results is 
worth the added computational effort only when the gust wave lengths of importance are 
comparable to the wing span or less. Although only the normal component of gust 
velocity is illustrated in Figure 1-1, lateral and longitudinal gust velocities are treated 
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Figure 1-1. Evolution of Gust 
Design Criteria 


t 


0(g) 

(t w 2 



WAVELENGTH A (FT.) 

J L L 


„-4 


-3 


„-2 




10"“* 10° 10" - 10' UT 

REDUCED FREQUENCY, g, RADIANS/FT. 

ISOTROPIC TURBULENCE; VON KARMAN 

<r 2. 1 + 8/3(1. 339LT2) 2 

<f>(Q) = r vi n/16 


[l + (1 . 3339LT2) 2 ] 


■ 0(g) 

<7 w 2 



I 


,-4 


.-3 


-L 


J 


„-2 


.-1 


10 “* 10 “ 10 “• 10 10 
REDUCED FREQUENCY, g, RADIANS/FT. 

ISOTROPIC TURBULENCE; DREYDEN 
a w 2 L 1 + 3L 2 g 2 


0(g) = 


(1 + g 2 L 2 ) 2 


Figure 1-2. Analytical Representation of 

Atmospheric Turbulence Spectra 

in the same manner. Longitudinal gust velocity effects on aircraft are usually considered 
negligible. • 

i 

Launch vehicle gust criteria development has paralleled the aircraft gust criteria in 
many ways. The initial emphasis was placed on determining the response to discrete 
gust velocities. The missile flight path is characterized by continuously changing 
velocity and altitude, whereas an aircraft is often at a fixed airspeed and altitude for 
long periods of time. The vertical path of. the missile through varying wind shears 
produces significant gust loads, whereas only the turbulence produced by wind shears 
affects the airplane loads. These considerations have led to the concept of computing 
gust loads on missiles due to flight through a wind-shear profile on a deterministic 
basis. While this method may have been adequate for missiles where fatigue life was 
not an important consideration, due to the large number of shuttle flights expected 
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it is suggested and recommended that gust loads be determined on a probabilistic 
basis for space shuttle. Although the method described would require that the flight 
path be segmented in sufficiently small increments to permit airspeed, altitude, and 
mass properties to be considered constant over the increment, this disadvantage ap- 
pears to be more than offset by the advantages of the probabilistic approach. In the 
present study, one-dimensional turbulence described by the Von Karman spectrum is 

assumed. 

The concept of computing response characteristics of the vehicle in the frequency 
domain is attractive from the unsteady aerodynamic complexity viewpoint. In general, 
it is easier to compute the unsteady aerodynamic pressure produced by oscillatory 
motion of the vehicle or oscillatory gust velocity than to compute the response of the 
vehicle to a step function change in downwash. (The latter step is the usual preliminary 
.step before applying Duhamel's integral for determining the response to an arbitrary 
variation in the downwash. ) By assuming harmonic motion, all recently developed 
oscillatory aerodynamic theories can be used to compute the frequency response func- 
tions. Subsequently, the frequency response functions can be combined with the random 
gust spectra to obtain response spectra, or they can be used to obtain a Fourier series 
solution of the response in the time domain to a discrete gust of arbitrary shape. Hence, 
solutions in both the frequency domain and time domain can be obtained. 
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SECTION 2 

TECHNICAL DISCUSSION 


This section presents a discussion of the equations of motion included in the computer 
program. The method is general, and may be applied to any vehicle whose motion is 
described in terms of modal generalized coordinates. 

2. 1 TURBULENCE RESPONSE 


Vehicle response is computed by solving sets of simultaneous equations obtained by 
applying Lagranges' equation of motion: 


d / d L \ 


dL SjQ - 

+ — Q 

dq £q r 

r r 


r = 1, 2, ... n 


where 

L = T - P 
T = kinetic energy 
P = potential energy 
# = dissipation function 
Q r = r^ generalized force 

XL 

q^ = r Ln generalized coordinate 
The deflection (h) at any point on the vehicle is given by: 


(2-1) 


n 

h (x,y,z,t) = y t h^ (x,y,z)q^(t) (2-2) 

r=l 

where h r is the normalized deflection at the point in the r^ mode shape. 

Using normal (with respect to the stiffness matrix) modes of vibration as generalized 
coordinates, the r^ 1 equation of motion is of the form: 


M q + D q + M 
rr r rr r rr 
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co q 
r r 



M 

rs 




2-1 


(2-3) 



where 

i M generalized mass 
rs 

D = generalized damping 
rr 

to = frequency of mode r 

6 * (0,r=s 

rs i l,r/s 


For the vehicle in the absence of a stability augmentation system, the generalized 
force includes the force produced by the gust and the force caused by the vehicle 
response. This can be written as: 

V 2X 5 rs +q Af <2 ' 4) 

s=l 

Response Gust 

where Q rg is the r^ 1 generalized aerodynamic force due to a unit value of the vehicle 
displacement in mode s, q^ is the gust amplitude, and Q r f is the r^ 1 generalized aero- 
dynamic force due to a unit gust. 


For harmonic motion, the equations of motion can be written in matrix form as: 


where 
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g r = structural damping coefficient for mode r 
co = gust frequency 

q r = amplitude of the r th generalized coordinate 
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(x,y,z)Ap s (x,y,z, co) dS 


( 2 - 8 ) 


where p is air density, b is the reference length, Ap is the net unsteady pressure 

O 

acting on the vehicle due to oscillation in mode s, and the integration is over the 
surface of the vehicle. Methods available for determining the unsteady pressure dis- 
tribution are described briefly in Section 2. 3. 


The generalized force due to the unit sinusoidal gust is given by: 


Q rf (<o) = 



(x,y,z)Ap f (x,y,z,co) dS 


(2-9) 


where Apj is the net pressure over the vehicle surface due to the downwash produced 
by a continuous sinusoidal gust wave traveling across the vehicle. This sinusoidal 
gust velocity may be expressed as: 


Wj = q^e 


icot 


= (cos oot + i sin cot) 

The downwash on an oscillating airfoil in a flow of velocity (V) is given by: 
w(t) = V a(t) - h(t) 


(2-10) 


( 2 - 11 ) 


where a is the angle of attack and h is the plunging velocity. For harmonic motion, 
the complex downwash may be written: 

w = V ca - i ooh (2-12) 


Equating 2-10 and 2-12 yields: 


W f = w 

- qf 

a = — cos co t 


- q f . „ 

h sin cot 

co 


(2-13) 
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If x* is the distance from a gust reference axis to the downwash at point i on the 
vehicle and k is the reduced frequency 



then Equations 2-13 and 2-14 become: 


- u - 

a. . = — cos k x. 
i V i 

q f b 

h i = " vk sin k x i 


(2-16) 

(2-17) 


where 



From the a and h distributions, the net pressure distribution (Ap^) due to the gust is 
obtained using the same unsteady aerodynamic theory employed to compute pressure 
distribution due to vehicle motion. 


Response transfer functions are readily obtained from the generalized coordinates 
computed in Equation 2-5. The acceleration transfer function is given by: 


H.(x,y,z, cu) 
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(2-18) 


Loading transfer functions can be obtained by modal displacement or load summation. 
The modal displacement technique is based on the concept that the deformed shape of 
a body defines the load (shear, moment, and torque) distributions. Assuming that n 
generalized coordinates are sufficient to define the deformed shape, the load transfer 
functions are of the form: 


1 t 

H(x,y,z,co) =— 23 F r (x,y,z) q r (w) (2 ' 1! 

qf r=l 

where F is the load (shear, moment, or torque) due to a unit value of. q r . 

r . 

The load summation method combines the loads due to inertia forces and air forces. 
The load transfer function is given by: 
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r=l 


F f (x,y,z, to) 

% 


( 2 - 20 ) 


th 


where F ' is the inertia load (shear, moment, or torque) due to the r degree of 
freedom oscillating at a unit frequency, F r is the load due to the aerodynamic pressure 
distribution developed by a unit amount of the r^ 1 generalized coordinate oscillating at 
frequency to, and F f is the load due to the aerodynamic pressure distribution produced 
by a gust wave of frequency to. 

2. 1. 1 RESPONSE TO RANDOM DISTURBANCE . The response power spectral density 
(PSD), 0 , is related to the gust PSD, 0 , and tye transfer function magnitude, | H | , 
by the equation: 


0 (co) = |H(oj)| 0 (co) 

X 


( 2 - 21 ) 


The two most common gust PSD forms are shown in Figure 1-2, where U - co/V, L is 
the scale of 
per second. 


the scale of turbulence measured in feet, and cr,„ is the rms turbulence velocity in feet 

W 


The ratio of rms response to rms gust velocity is given by: 

1/2 


A = 


co. 


/ 


0 (co) dto 
x 


( 2 - 22 ) 


where oo n is the integration cutoff frequency. The "characteristic frequency" of the 

V 

response, or number of zero crossings with positive slope per second, is expressed 
by: j 


N = 


0 27TA 


; co. 


/ 


co 0 (co) dto 

X 


L o 


1/2 


(2-23) 


The total number of response peaks per second is given by: 

ll/2 


Nl AN 


co. 


r c 4 

/ co 0„(co) dto 

Lo 


x 


(2-24) 
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9.. 1.2 RESPONSE TO DISCRETE DISTURBANCE. The Fourier series technique is 
used to "Resent the discrete forcing function. This is written: 


M 


f(t) = a o + S (a m C0S ^m 1 + b m sin 
m=l 


(2-25) 


where 


f(t) is the forcing function 

M is the number of terms in the series 

t = time 

co = m tt/T 
m 


i 

a = — f f(t)dt 

a o 2T 1 


m 


= — f f(t) cos to t dt 


-T 


m 


4 / 

-T 


f(t) sin w m t dt 


The resulting force time-history f(t) is periodic with period 2T. 


As 

to 


prescribed by Equation 2-25, the force is stationary. It is often more 
think of the vehicle as being fixed, with the disturbance moving over it. This is 


111 liUlUV y ^ . • 1J 

described by substituting (x/V-t) for t in Equation 2-25, which yields: 


M 


f(t) = a + 53 a m cos (k m x- w m t) + b m sln(k m _x- w m t) 


(2-26) 


m=l 


k = cu m b/V 


m 


x = x/b 

b is the reference length 


where 



V is the vehicle velocity 

x is the distance from a point on the vehicle to the gust reference axis 


By considering that the cosine terms are symmetric and the sine terms are antisym- 
metric, Equation 2-26 can be written: 

M 

£(t) = a o + X) f a m cos( “m t - k m 5) ‘ b m sln <“m t - k m : >l (2 ‘ 27) 

m=l L J 

In complex form, Equation 2-27 can be written as: 


f(t) = a + Re 
' o 

1 

i(a! m t-k x) 

a e 
m 

-Im 

‘ M i(w m t-k x) 

L b m e 


m=l 



m=l 


The vehicle response to this disturbance is given by 


Z W= Z r( a m H R m + b m H I m ) cos ^rn t + ( b m H R m ~ a m H Im ) sin w m t] (2-29) 
m=l L J 

where H R and H T are the real and imaginary components of the m™ transfer function 


v m a m 
evaluated at frequency co m , i. e. : 


H = H p + iHr 
m »m L 


m 


(2-30) 


2. 2 STABILITY AUGMENTATION SYSTEM SIMULATION 


For a vehicle having a stability augmentation system (SAS), Equation 2-5 is modified 
to include the control surface aerodynamic and inertia forces so that the equation 
becomes: 


[A(co)] jq(w)j = - }A f (co)J - Mg + Q § (co)J 6(co) 


(2-31) 


th 


where M r . is the generalized mass term coupling the r vehicle mode to the control 
o 

surface unit rotation, and Q r ^ is defined by: 

/h r (x,y, Z )Ap 6 (x,y,z, W )dS 

4pb. w s 
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(2-32) 





where the integration is over the entire vehicle surface and Ap 6 is the net unsteady 
pressure distribution acting on the vehicle due to unit amplitude oscillation of t e 

control surface at frequency co. 

It Is interesting to note that thevehicle response to a time-dependent control surface 
motion can be obtained by applying the Fourier series method outlined in Section 2. 1. 2 
to Equation 2-31 to describe 6(t) and setting Af equal to zero. 

The general block diagram for a vehicle with a one-channel SAS is shown in Figure 2-1. 



Figure 2-1. SAS Block Diagram 


Control surface rotation is given by: 

6(a)) = -R(uj) - Hp(a>) H-^(to) Z^(co) 


(2-33) 


where R(u)) is the control input, H F and H 1 are system transfer functions, 
the response. The response is given by: 


End is 


Z 1 ( 01 )= 52 h r (1) q r (w) (2 34 ? 

r=l 

where hjP 0 is the normalized deflection at point (1) in the r & mode shape and 5 r (co) is 
the r th generalized coordinate evaluated at frequency co. Substituting Equation 2-34 
into 2-33 yields in matrix form: 


6 (co) = -R(co) - [g(co)J jq(co) 

where 

|G(m)j =H F (U)H l (m)jh< 1 )J 


(2-35) 


(2-36) 
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Substituting Equation 2-35 into 2-31 yields the following expression for the generalized 
coordinate response. 


;j = -[[ A] - jp 6 J lgj] 


-1 


where 


F 6 (co)| = [m 6+ Q 6 (o.)| 


Afj - | F 6| R i 


(2-37) 


(2-38) 


It can be seen from Equation 2-37 that the effect of an automatic control system is to 
modify the A- matrix. For a vehicle with an SAS and no control input command (R = 0) , 
Equation 2-31 may be written: 


-Krh 


(2-39) 


where 


[a*] - [A] - 


F 6 j LGj 


Equation 2-39 has exactly the same form as Equation 2-5, the response of the unaug- 
mented vehicle. 

For a system with multiple sensors such as shown in Figure 2-2 , the expression for 
g(co ) becomes: 

m 

[G(w)J =h f m E (2-40) 
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For a system having m independent control surfaces, the generalized coordinate 
response is given by: 



In other words, additional control channels are additive to the single control surface 
case given by Equation 2-37. 

2.3 AERODYNAMICS 

Parallel-mounted vehicles, such as those being proposed for the space shuttle, present 
unique problems to the aeroelastician. Interference effects between the booster and 
orbiter vehicles in the mated configuration have been found in wind tunnel tests to be 
significant, particularly in the yaw plane. Figure 2-3 shows the derivatives of the 




Figure 2-3. C n ^ and Cy^ Versus Mach Number for Booster Showing 

Effects of Interference Effect Due to Orbiter 
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yawing moment (C n ^) and the side force coefficient (C y ^) as a function of Mach number 

for the configuration shown in’ Figure 2-4. Interference effects tend to double the 
isolated values at the lower Mach numbers. 


The space shuttle is characterized by a large body relative to the wing and tail surfaces. 
This is because lifting surfaces are required only when the propellant is depleted and 
the vehicle is light. Figure 2-5 shows the derivative of normal force coefficient (C^) 

as a function of Mach number for the booster. The top curve is a plot of wind tunnel 
data for the complete booster, the middle curve for the same configuration with the 
wings removed, and the bottom curve the difference between the two. Subsonically, 

Cn for the body-tail is about 85 percent of the total. Transonically, the body-tail 

is about half of the total, increasing to about 70 percent of the total supersonically. 

The nature of the space shuttle mission is such that it flies at subsonic, transonic, 
supersonic, and hypersonic speeds. Figure 2-6 shows a typical plot of dynamic pres- 
sure versus Mach number during boost and entry flight. Angle of attack during boost 
flight ranges from about plus 5 degrees to minus 10 degrees. The separated vehicles 
enter at very high angle of attack and then transition to airplane-type flight at low 
supersonic or subsonic speeds. The requirement for subsonic cruise to the landing 
site is the reason that the airfoils tend to be thick (10 to 20 percent) and have subsonic 
planforms. 



Figure 2-4. Typical Space Shuttle Configuration 



DYNAMIC PRESSURE (psf) 




Figure 2-6. Typical Space Shuttle Booster Dynamic Pressure Versus Mach Number 
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From the preceeding discussion, it is obvious that any aerodynamic theory employed 
to predict space shuttle response to atmospheric disturbance should account for body 
aerodynamics, aerodynamic interference effects, thickness effects, and angle of attack. 

A survey of available unsteady and quasi-steady aerodynamic theories has been conducted 
to determine their applicability to the space shuttle aeroelasticity problem. Theories 
investigated include: 

a. Unsteady theories 

Subsonic 

Strip theory (Reference 2) 

Kernel function (Reference 3) 

Doublet lattice (Reference 4) 

Transonic 

Transonic box (Reference 5) 

Supersonic 

Mach box (Reference 6) 

Hypersonic 

Piston theory (Reference 7) 

b. Quasi-steady methods 

Woodward method (Reference 8) 

Modified strip theory (Reference 9) 

Table 2-1 summarizes the applicability of these theories when the following effects 
are significant. 

a. Large body with small aerodynamic surfaces 

b. Airfoil thickness 

c. Body/body and body/wing interference 

d. Control surfaces 

e. Angle-of-attack 

f. High (low) aspect ratios 

Table 2-1 indicates that the present state-of-the-art in unsteady aerodynamics is in- 
adequate for analyzing complex configurations such as space shuttle. The Air Force 
and NASA are continuing to sponsor research in this area and hopefully the picture 
will brighten in the future. For the meantime, however, when interference effects, 
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Table 2-1. Applicability of Referenced Theories to Significant Effects 



Theory 

a 

b 

Effect 
c d 

e 

f 

Subsonic 

Strip Theory 

No 

No 

No 

Yes 

No 

M-H 


Kernel Function 

No 

No 

No 

Yes 

No 

All 


Doublet Lattice 

No 

No 

No 

Yes 

No 

All 

Transonic 

Box 

No 

No 

No 

Yes 

No 

L-M 

Supersonic 

Mach Box 

No 

No 

No 

Yes 

No 

L-M 

Hypersonic 

Piston Theory 

Yes 

Yes 

No* 

Yes 

Yes 

All 

Quasi-steady 

Woodward 

Yes 

Yes 

Yes 

Yes 

Yes 

All 


Modified Strip Theory 

No 

No 

No 

Yes 

No 

M-H 


* Interference effects become small hypersonically 


thickness, angle -of-attack, etc. are significant factors and the reduced frequency 
range of interest is not so high that unsteady lag effects dominate, it is recommended 
that Woodward's steady aerodynamic method be employed at speeds below hypersonic. 
Piston theory appears adequate in the hypersonic Mach number range. 

The common approach when using steady-state aerodynamics to derive generalized 
forces for flutter and dynamic response problems is to compute the complex downwash 
at a point j due to a unit amount of the k th generalized coordinate from the equation: 




+1 “k h jk 


(2-42) 


where the first term is the downwash due to the flight velocity, V, and the modal angle 
of attack, a, and the second term is the downwash induced by the surface plunging 
motion of amplitude h and frequency to. 

The steady-state method of Woodward was applied by Brignac and Shelton (Reference 
10) to flutter analysis of the F-lll horizontal stabilizer. Analytical and experimental 
results were compared for flutter Mach numbers of 1. 09 and 1.37. The results are 
shown in Tables 2-2 and 2-3. At Mach 1. 09, where no applicable unsteady aerodynamic 
theory was available, the predicted flutter velocity using Woodward aerodynamics was 
8 percent low and the flutter frequency was 2 percent high. This is remarkable corre- 
lation considering the transonic Mach number and the relatively high reduced frequency 
at flutter (about 0. 8 based on the mean semi-chord). At Mach 1.37, the predicted 
flutter speed using the Woodward method agreed exactly with test results. The pre- 
dicted flutter frequency was 16 percent high. Unsteady Mach box theory applied to 
this condition yielded a 3 percent higher flutter speed. 



J 


Table 2-2. Comparison of Quasi-Steady and Experimental 
Flutter Results at Mach 1. 09 



Mach 

No. 

Mass -Density 
Ratio (ji) 

V F 

Vp (test) 

f F 

f-p (test) 

Test 

1.09 

13. 10 

1.00 

1. 00 

Analysis 

(Quasi-Steady) 

1.09 

13. 10 

0.92 

1. 02 

Table 2-3. 

Comparison of Quasi-Steady, Unsteady, and 
Experimental Flutter Results at Mach 1. 37 


• 

Mach 

Mass-Density 

V F 

f F 


No. 

Ratio (\j) 

Vp (test) 

fp (test) . 

Test 

1.37 

12.88 

1.0 

1. 00 

Unsteady 
(Mach Box) 

1.37 

12. 88 

1.03 

1.00 

Quasi-Steady 

(Woodward) 

1.37 

12.8 

1.00 

1. 16 


It becomes evident that in some cases the use of steady-state aerodynamics yields 
more accurate flutter calculations than does unsteady aerodynamic theory. Apparently, 
unsteady lag effects are sometimes less important than the proper definition of steady- 
state effects due to interference, thickness, angle of attack, etc. The problem then 
becomes one of defining the steady-state aerodynamic pressure distribution. As shown 
in Figures 2-7 through 2-9 taken from Reference 11, the Woodward method appears to 
have excellent potential in this regard. 

2.4 PROPELLANT SLOSHING AND ENGINE DYNAMICS 


Propellant sloshing is accounted for using the method presented in Reference 12, which 
defines a "sloshing mass" distribution along the length of the tank. A spring-mass 
representation is applied at the center of gravity of the sloshing mass. Spring stiffness, 
k g , is given by: 


k s = m s co s 


(2-43) 


where m g is the sloshing mass and to g is the slosh frequency. 
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Fuselage- Fin Configuration (a = 4 deg, Mach 0. 8) 



Fin Configuration (a = 4 deg, Mach 0. 8, $ = 0 deg) 



Figure 2-9, Fin Pressure Distribution for B-58 Wing- Fuselage- 
Fin Configuration (a = 4 deg, Mach 0. 8, /3 = 0 deg) 
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There are two methods for including sloshing in the vehicle dynamics. The slosh spring- 
mass system can be included in the mathematical model used for the vehicle modal an- 
alysis, or it can be mass-coupled to the normal modes derived based on rigid propellant. 
In the latter technique, the terms in the generalized mass matrix that couple the vehicle 
modes to the slosh modes are given by: 

M = m 0 . ( 2-44) 

is s is 

where M^ s is the generalized mass term coupling the ft 1 vehicle mode to the slosh mode 
and 0-[ s is the deflection of the slosh mass due to vehicle mode i. 

The thrust force acting normal to the booster centerline (assuming small angles) is: 

F t = -T (0 +6) (2-45) 

where T is the engine thrust, 0 is the slope at the gimbal point, and 6 is the engine 
gimbal angle relative to the deformed centerline as shown in Figure 2-10. 



The slope of the deformed centerline at the gimbal point is given by: 


n 


8(t) = £ h r ' 'J r < t) 

r=l 


(2-46) 


M ttl 

where h ' is the r tn modal slope at the gimbal point and is the r generalized co- 
ordinate. | 

Substituting Equation 2-46 into 2-45 yields: 


n 


F T = - T (£ h r'\ + 6 


,r=l 


(2-47) 


The i^ 1 generalized force due to thrust is given by: 
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Q Ti = - T (S h r /q r + 6 ) h i 


(2-48) 


»>r=l 


where ^ is the 1 th modal deflection at the gimbal point. • 
The inertia force at the engine eg due to engine gimballing is: 


F = -m 4 6 

where m is the mass of the engine 
engine eg. 


(2-49) 


and & is the distance from the gimbal point to the 


The inertia moment about the engine eg is: 

Mj. =16 

where I is the engine mass moment of inertia about the engine eg. The generalized 
inertia force acting on mode i is written: 


(2-50) 


V f i(V*V) +m iV 

Substituting Equations 2-49 and 2-50 into 2-51 yields: 

Qj. = Ijm-^h. + + l)h/ J 6 

which can be rewritten: 

Q Ii = (-m^h.+Igh. / )*6 

i 

where I is the engine mass moment of inertia about the gimbal point, x. e. : 

s 

I = m £ 2 + I 
g 

The total generalized force acting on mode i due to the gimballed engine is then: 

/ n \ 


(2-51) 


(2-52) 


(2-53) 


Q. = Q Ti + %-~ T (S . .V \ + 6 j h i 


+ (-m^hi + Igh^) *6* 


(2-54) 
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Assuming harmonic motion 


and dividing Equation 2-54 by 4 p b 3 co 2 yields: 


Q.(co) 


Q^(co) 
4pb 3 co 2 


1 

4p b 3 


n 


T 23 VM^i 

+ (mJh 1 -I g h i / 


co 



(2-55) 


which can be written in matrix form as : 



6 x (co) + 




where: 



(2-56) 


6 l( co)= Lh'jjq(co)J 

4pb co 2 

4 p b v w ' 

6 2 (co) = 6 (co) 

\ 

Incorporating these expressions into the equations of motion parallels exactly the 
method outlined in Section 2.2, 



SECTION 3 


COMPUTER PROGRAM DESCRIPTION 


A computer program, based on the equations presented in Chapter 2, has been written 
to calculate the response of an elastic aerospace vehicle to random and discrete turbu- 
lence. The program, documented in Reference 1, is written in Fortran IV for the 
UNIVAC 1108 and CDC 6400 computers. 

Inputs to the program consist of modal data, generalized aerodynamic forces, control 
system description, flight data, and discription of the disturbing function. Program 
outputs include tabulations of A and N Q and SC 4020 plots of transfer functions, response 
power spectral densities, and Nyquist stability traces. Response time histories are 
plotted for the discrete turbulence case. The program is dimensioned to have the 
following capabilities. 

20 generalized coordinates, 

20 structural response parameters, 

20 reduced frequencies, 

10 interpolated points between input reduced frequencies 
(for a total of 210 response frequencies). 

The following features add to the usefulness of the program. 

a. Most large data blocks can be input or output on magnetic tape. 

b. Most of the equations contain input scalar multipliers. 

c. Several standard input power spectral density functions are available in the program; 
in addition, the user may construct his own using straight-line segments. 

d. Several common time-varying forces are available or the user may construct his 

own. i 

e. Extensive use is made of the SC 4020 plotter. 

f. A comprehensive error routine assists in finding input errors^ 

g. Provisions are made to compute quasi-steady forces from steady-state aerodymamic 
input data. 

h. Ten types of control system elements are available for defining the stability aug- 
mentation system (SAS). 



The digital computer program, is divided into six main parts: 

Part 0 computes the dynamic matrix representing the vehicle SAS. 

Part 1 interpolates the aerodynamic terms and calculates the response of the 
generalized coordinates. 

Part 2 computes the transfer functions. 

Part 3 determines the response to a random input. 

Part 4 calculates the response to a discrete input. 

Part 5 plots the results of Parts 2, 3, and/or 4 on the SC 4020 plotter. 

The program is set up such that Part 0 is run independently and Parts 1 through 5 are 
run independently or in series. 

The input to the program can consist of problem cards and magnetic tape. Program 
output consists of user-specified combinations of printed information, punched cards, 
magnetic tape, and SC 4020 plots. 

Section 4 presents the results of an application of this program to the turbulence re- 
sponse analysis of a typical space shuttle configuration. 



SECTION 4 


PRELIMINARY RESULTS OF SPACE SHUTTLE 
TURBULENCE RESPONSE ANALYSIS 


A typical space shuttle configuration was analyzed to determine its response to random 
turbulence. The configuration consists of the delta wing booster and the straight wing 
orbiter shown in Figure 4-1. Dui’ing ascent flight> prior to separation, the orbitei is 
carried "piggyback" by the booster. 

Symmetric and antisymmetric responses were computed at a point on the ascent tra- 
• jectoi’y corresponding to 0. 9 Mach and 22, 800 feet altitude. Symmetric response of the 
booster during flyback was determined at 0. 9 Mach and 40,000 feet altitude. 


For both the ascent and the booster return flight analyses, the Von Karman gust spec- 
trum was assumed with a 2500-foot scale of turbulence. 

The analyses were conducted with the stability augmentation system (SAS) both opera- 
tive and inoperative. Autopilot gains and sensor locations were assumed without the 
aid of an elastic vehicle control stability analysis. The unaugmented vehicle results 
are therefore more realistic. No attempt was made in the present study to optimize 
the -SAS. In some SAS-on cases, instabilities are suggested by plots of the stability 
determinant. 


Vibration modes were computed using a finite element method developed at Convair 
Aerospace. For these mode shapes, unsteady aerodynamic forces were generated 
using the kernel function theory (Reference 3) for the aerodynamic surfaces and 
slender body quasi-steady aerodynamics (Reference 7) for the body. 


4. 1 BOOSTER RETURN FLIGHT 


Referring to Figure 4-2, the rigid body equations of motion are as follows for transient 
pitch plane disturbances: 


Z - 0V = - 


QSg 

W 


( c I tt “ + 0 L 6E « E ) 


(4-1) 


0 = 


QSc 


( C Ma a 


-M^e) 


+ Cm 6i 


(4-2) 


,4-1 


I 







e 



where 


Figure 4-2. Pitch Plane Sign Convention 

Z is normal acceleration of the eg (ft/sec 2 ) 

"(5 is pitch acceleration (rad/sec 2 ) 

Q is dynamic pressure (lb/ft 2 ) 

S is reference area (ft 2 ) 
g is acceleration due to gravity (ft/sec 2 ) 

W is weight (lb) 
c is reference length (ft) 

I is mass moment of inertia about the eg (slug-ft 2 ) 

Cl is the derivative of the lift coefficient with respect to angle of attack 
(per rad) 

Cl^ is the derivative of the lift coefficient with respect to elevon deflec- 
E tion angle (per rad) 

Cmq, * s derivative of the pitching moment coefficient with respect to 
angle of attack (per rad) 

C M; . is the derivative of the pitching moment coefficient with respect to 
elevon deflection angle (per rad) 

V is the flight velocity (ft/sec) 

is elevon deflection angle (rad) 

a is the angle of attack (rad) 


Z 

a = — + a , 
V gust 


4-3 



J 

The booster flyback pitch plane SAS block diagram is shown in Figure 4-3. The con- 
trol equation is : 

T A®E + 6 E =K A <K 8 9 + K e®" K Z 2) ' <4 " 3> 


where 


K 9 

K e 


is the elevon actuator time constant (sec) 

is the elevon actuator gain (rad /rad) 

is the pitch displacement gain (rad/rad) 

is the pitch rate gain (rad/rad/sec) 

is the normal acceleration feedback gain (rad/ft/sec ) 


It is assumed that the canard is used for trim only and that the SAS sensors are located 
at the eg. 



Figure 4-3. Booster Flyback Pitch Plane SAS 


Numerical values for the booster flyback configuration are listed in Table 4-1. 

The booster finite element model used for the modal analysis represents the aerodynamic 
surfaces as a redundant network of spars and ribs and the fuselage as a 12 grid-point 
beam. Canard idealization is shown in Figure 4-4. Beam segments represent spars 
and ribs. Plate elements are used for the cover skins. Solid lines represent locations 
of physical members and dashed lines indicate members added for analysis purposes. 

The solid circles show locations of masses in the simulation. The member extending 
inboard from grid point 167 is the pivot shaft. 

; • 4-4 







Table 4-1. -Data for the Booster Return Flight Analysis 


Mach Number 

0.9 

% 

2. 3 per rad 

Altitude 

40,000 ft 

■°“o 

-0. 029 per rad 

Q 

224 lb/ft 2 

Cl «e 

0. 435 per rad 

V 

874 ft/sec 

Cm 6. E 

-0. 326 per rad 

s 

7791 ft 2 

t a 

0. 0667 see 

c 

71.69 ft 

Ka 

4. 0 rad/rad 

w 

694,435 lb 

K e 

1. 0 rad/rad 

I 

84. 307 X lO 6 slug-ft 2 

K 8 

1. 0 rad/rad/sec 

eg Location 

2013 inches aft of nose 

K z 

0. 124 rad/ft/sec 2 


Wing idealization is shown in Figure 4-5. Grid points 124, 127, 130, 133, and 136 
represent the airbreathing engines. The actual wing is covered with corrugated skis. 

It was assumed that the skin was fully ineffective in bending and only partially effective 
in carrying shear. The idealized wing is suspended from the fuselage by springs located 
at grid points IS, 21, 23, and 27. 

The first eight symmetric mode shapes for the booster return configuration are shown 
in Figures 4-6 through 4-13. The first mode is wing bending, and the second and third 
are combinations of wing torsion and fuselage bending. The fourth mode is canard pitch. 

Results or the turbulence response analysis for booster return symmetric flight without 
the SAS are summarized in Table 4-2. For the one-degree-of-freedom (DOF) analysis, 
only rigid-body vertical translation is considered. For the two-DOF analysis, the rigid 
body pitch mode is included. The 10-DOF case considers the two rigid body modes and 
first eight elastic modes of the vehicle. As can be seen in the table, elastic modes are 
required to predict the response properly^ 

» 

Figure 4-14 shows A and N Q for normal acceleration at points along the fiiselage. 

Except for the nose region, the two-DOF case adequately predicts the A distribution 
but not the N 0 values. 

To lend more physical insight to the results shown in Figure 4-14, assume a design 
rms gust velocity of 50 feet per second. The rms normal acceleration at the booster 
nose is then 0. 5 g. 

Acceleration transfer functions for three points on the fuselage are shown in Figures 
4-15 through 4-17. Corresponding PSD plots are given in Figures 4-18 through 4-20. 
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Figure 4-4. Booster Canard Structural Idealization 
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Table 4-2. Booster Return Flight Load Summary (SAS Off, Symmetric) 



1 DOF 

2 DOF 

10 DOF 

Response Item 

A 

N 0 (Hz) 

A 

N 0 (Hz) 

A 

N q (Hz) 

Wing Root Shear (lb) 

1008. 0 

0.780 

1331. 1 

0. 577 

1174. 6 

1.217 

Wing Root Bending 
Moment (in-lb) 

290,000 

0.835 

382,810 

0. 622 

334,730 

1.318 

Wing Root Torque About 
FS 3488. 5 (in-lb) 

41,424 

2.789 

48,819 

2.38 

157. 760 

3.828 

Fuselage Shear at 
FS 2800 (lb) 

729. 6 

0.983 

915.28 

0.964 

843.5 

1. 422 

Fuselage Bending Moment 
at FS 2800 (in-lb) 

350,000 

1.511 

251,780 

1. 181 

413,210 

3.829 


The transfer functions for fuselage shear and bending moment at fuselage station (FS ) 
2800 are shown in Figures 4-21 and 4-22, respectively. The corresponding PSD 
curves are contained in Figures 4—23 and 4-24. Wing root shear, bending moment, 
and torque about FS 3488. 5 are shown in Figures 4-25 through 4-27, respectively. 
Figures 4-28 through 4-30 are the corresponding response PSD plots. 

An indication of the vehicle stability with SAS off is shown in Figures 4-31 and 4-32. 
Figure 4-31 is a plot of the magnitude and phase angle of the determinant of the A- 
matrix evaluated over the frequency range of interest, where the A-matrix is defined 
by Equation 2-5a. A polar plot of the same information is given in Figure 4-32, which 
is analogous to a Nyquist stability diagram. The continuous counterclockwise encircle- 
ments of the origin indicate that the system is stable throughout the frequency range 
(Reference 13). 

Calculations with the SAS on considered the rigid body vertical translation and pitch 
modes and the first six elastic modes. Acceleration transfer functions for the three 
fuselage points are shown in Figures 4-33 through 4-35. The corresponding PSD 
curves are given in Figure 4-36 through 4-38. Figure 4-39 shows the A and N 0 values 
for fuselage normal accleratiori versus fuselage station. Comparing the accelerations 
with those shown in Figure 4-14 for SAS off indicates that SAS suppresses the elastic 
mode response. With SAS on, nose acceleration is reduced to one-half the SAS-off value 
and the characteristic frequency over the forward three-fourths of the vehicle is re- 
duced considerably. 

The SAS-on A-matrix determinant characteristics are shown in Figure 4-40 and 4-41. 
Erratic behavior in the range of 4 to 4. 5 Hz suggests a system instability. 
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Frequency = 5. 47 Hz 



Figure 4-8. Booster Return Flight, Third Symmetric Mode 


Frequency = 6.28 Hz 



Figure 4-9. Booster Return Flight, Fourth Symmetric Mode 
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Figure 4-15. Booster Return Flight Fuselage Acceleration Transfer 
° Function (SAS Off, Symmetric) at Crew Station 


PHASE ANGLE ACCELERATION CG) PER 1.0 FT/SEC GUST. 6/FPS 
(DEGREES) 




ACCELERATION C6) PER 1.0 FT/SEC GUST» 6/FPS 



Figure 4-17. Booster Return Flight Fuselage Acceleration Transfer 
Function (SAS Off, Symmetric) at Lfhj Tank 
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Figure 4-18. Booster Return Flight Fuselage Acceleration 
PSD (SAS Off, Symmetric) at Crew Station 
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Figure 4-19. Booster Return. Flight Fuselage Acceleration 
PSD (SAS Off, Symmetric) at LOg Tank 
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Figure 4-20. Booster Return Flight Fuselage Acceleration 
PSD (SAS Off, Symmetric) at LHg Tank 
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Figure 4-22. Booster Return Flight Fuselage Bending Moment 

Transfer Function at FS 2800 (SAS Off, Symmetric) 
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Figure 4-26. Booster Return Flight Wing Root Bending Moment 
Transfer Function (SAS Off, Symmetric) 
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Figure 4-27. Booster Return Flight Wing Root Torque About 

FS 3488. 5 Transfer Function (SAS Off, Symmetric) 
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Figure 4-29. 


Booster Return Flight Wing Root Bending 
Moment PSD (SA S Off, Symmetric) 
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Figure 4-31. Booster Return Flight A-Matrix Determinant Magnitude 
and Phase Angle (SAS Off, Symmetric) 
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Figure 4-32. Booster Return Flight A-Matrix Determinant 
Polar Plot (SAS Off, Symmetric) 
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Figure 4-33. Booster Return Flight Fuselage Acceleration Transfer 
Function (SAS On, Symmetric) at Crew Station 
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Figure 4-34. Booster Return Flight Fuselage Acceleration Transfer 
Function (SAS On, Symmetric) at LOg Tank 
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Figure 4-35. Booster Return Flight Fuselage Acceleration Transfer 
Function (SAS On, Symmetric) at LH 2 Tank 
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Figure 4-36. Booster Return Flight Fuselage Acceleration 
PSD (SAS On, Symmetric) at Crew Station 
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Figure 4-37. Booster Return Flight Fuselage Acceleration 
PSD (SAS On, Symmetric) at LO£ Tank 
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Figure 4-38. Booster Return Flight Fuselage Acceleration 
PSD (SAS On, Symmetric) at LH^ Tank 
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Figure 4-41. Booster Return Flight A-Matrix Determinant 
Polar Plot (SAS On, Symmetric) 
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4. 2 ASCENT FLIGHT 

1 

Referring to Figure 4-42, the rigid body symmetric equations of motion are: 


Sum of forces in z-direction: 

SF, - - C La « QS + T 6 pitch . (4-4) 

Sum of moments about y-axis: 

SMy = Cj^aQS (5 cg - ; AC ) + T 6 pltch (Sq - x cg ) , (4-5) 

Pitching acceleration (about y-axis): 

q^EMy/Iyy (4-6) 

Acceleration in z-direction: 

w = 2F Z /M + qu (4-7) 


In keeping with standard aircraft notation, u, v, and w are the velocity components in 
the x, y, and z directions respectively and p, q, and r are the three angular velocity 
components. 



Figure 4-42. Ascent Symmetric Flight Sign Convention 


The pitch plane SAS block diagram is shown in Figure 4-43. 

The rigid body equations of motion for transient antisymmetric flight using the sign 
conventions in Figure 4-44 are: 

Sum of forces in y-direction: 

SF y = Cy^QS- T 6 yaw , (4-8) 

4-41. 




Figure 4-43. Ascent Symmetric Flight SAS Block Diagram 




Figure 4-44. Ascent Antisymmetric Flight Sign Convention 
Sum of moments about x-axis: 

2“x = [( % + Ci 6a 6 a ] QS^ roU T6 roll 

Sum of moments about z-axis: 

£m z = (c n(J B ref + Xjj ^QS + T 6 y aw (; G - J cg ) 
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(4-9) 


(4-10) 





(4-11) 


Rolling acceleration (about x-axis): 



Yawing acceleration (about z-axis): 

r = (£M +1 p\/I ' ( 4 “ 12 ) 

z xz / zz 

Acceleration in y-direction: 

v = /M - ru ( 4 “ 13 ) 

y 

In Equation 4-9, 6 a is the aileron deflection defined such that a positive value causes 
positive roll. 

The SAS block diagram for ascent antisymmetric flight is shown in Figure 4-45. 



Figure 4-45. Ascent Antisymmetric Flight SAS Block Diagram 
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Engine thrust and inertia force equations are derived in Section 2. 4 for bending modes. 
The generalized forces produced by engine roll torques are given by: 

Q=X' 0j(m^6 -T5 n \ 

i roll x i \ roll roll/ 

where £ ro ll is the effective roll thrust moment arm and 0 X . is the slope about the x 
axis at the gimbal point due to mode i. 

The propellant sloshing analysis results for 59 seconds flight time are summarized in 
Table 4-3. 



Table 4- 

-3. Sloshing Parameters 


Tank 

Slosh Mass 
(lb-sec 2 /in.) 

Slosh Spring 
(lb/in.) 

Slosh eg (Booster Station) 
(in.) 

lh 2 

63.7 

373 


2863 

lo 2 

1030 

6020 


1622 


Basic data for the symmetric and antisymmetric ascent analyses are listed in Tables 
4-4 and 4-5, respectively. 


Table 4-4. Basic Data - Symmetric Ascent Flight 


t = 59 sec 

0 

T = 6. 3 x 10 lb 

Mach number = 0. 9 

Q = 486 lb/ft 2 

Altitude = 40,000 ft 

S = 7797 ft 2 

M = 1. 112 X 10 5 slugs 

Cl = 3. 84 per radian 

1^ = 411 X 10 6 slug-ft 2 

K a = 1 rad/rad 

x = 183. 1 ft 

T A = 0. 0667 sec 

X AC “ 213. ^ ft 

Kg = 1. 0 rad/rad/sec 

x Q = 312.5 ft 

Kg = 1.4 rad/rad 


For each of the 12 gimballed engines, the mass properties are: 
Mass = 165 slugs 

2 

Mass moment of inertia about gimbal point = 2505 slug-ft . 
Center of gravity distance aft of gimbal point = 20. 5 inches 
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Table 4-5. Basic Data - Antisymmetric Ascent Flight 


t = 59 sec 

Mach number = 0. 9 

Altitude = 40,000 ft 

5 

M = 1. 112 X 10 slugs 
I = 32.9 X 10 6 slug-ft 2 

XX 

T zz = 382 x 106 slu g _ft2 

I = -13. 7 X 10 6 slug-ft 2 
xz 

X cg = 183.1 ft 

x Q = 312.5 ft 

T = 6. 3 X 10 6 lb 

Q = 486 lb/ft 2 

S = 7797 ft 2 

j l . = 15 ft 
roll 

C v = -2. 40 per radian 
y P 


C 1 = -0. 424 per radian 

P 

C j> . = 0. 0658 per radian 

6 a 

C n ^ = 1. 32 per radian 

B ref - 145 ft - 

Z 1 = -6.34 ft 

X x = 99. 7 ft 

K 1 = 1 rad/rad 

T 1 = 0. 0667 sec 

Kg = 3. 0 rad/rad 

Tg = 0. 0667 sec 

K . = 1. 0 rad/rad 

0 

=0.7 rad/rad/sec 

K , = 1. 4 rad/rad 

¥ 

K; = 1.0 rad/rad/sec 

<P 


The booster structural idealization for the modal analysis was the same as that for the 
booster return flight symmetric analysis, with the addition of the vertical tail for the 
antisymmetric analysis. Tail idealization is shown in Figure 4-46, where the solid 
circles represent locations of masses. 

The rigid orbiter was attached to the booster by springs located at Stations 1829 and 
2657, as indicated in Figure 4-47. It was assumed that all axial load is transmitted at 
the forward attachment point, as is all yawing moment. 

The first six symmetric modes for ascent flight at 59 seconds are shown in Figures 
4-48 through 4-53. The orbiter center line is depicted by the line parallel to and above 
the booster fuselage center line. The first mode is body first bending; the second is 
wing first bending. The orbiter pitch mode is the fourth symmetric mode. 

Figures 4-54 through 4-60 show the first seven antisymmetric modes. The fundamental 
mode is orbiter yaw, and the second and third are first body bending and first wing 
bending, respectively. The fourth mode is tail first bending. 
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Figure 4 - 46 . Booster Vertical Stabilizer Structural Idealization 
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TYPICAL CROSS-SECTION LOOKING AFT 


The following stiffnesses relate Point 2 to Point 3. 



FWD (STATION 1829) 

AFT (STATION 2657) 

K x 

2.06X10 6 

0 lb/in. 

Ky 

0.223 X 10 6 

0. 274 x 10 6 lb/in. 

K z 

2.65 X 10 6 

1. 06 x 10 6 lb/in. 

K 0x 

5.97 X 10 6 

2.56 x 10 10 in.-lb/rad 

K e, 

0 1 

0 in. -lb/rad 

y 

K0 z 

4. 50 X 10 10 

0 in. -lb/rad 

Note: Links 

j : . 

1-2 and 3-4 are assumed rigid 



Figure 4-47. Booster /Orbiter Attachment Stiffnesses 
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Frequency = 6. 18 Hz 



Figure 4-52. Ascent Flight Fifth Symmetric Mode 



Figure 4-53. Ascent Flight Sixth Symmetric Mode 
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Figure 4-55. Ascent Flight Second Antisymmetric Mode 
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Figure 4-58. Ascent Flight Fifth Antisymmetric Mode 
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Figure 4-60. Ascent Flight Seventh Antisymmetric Mode 
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Symmetric turbulence response analysis results tor the ascent configuration at t = 59 
seconds are summarized in Table 4-6 tor the SAS-off condition. The 3-DOF case 
considers rigid body vertical translation plus the two slosh modes. Rigid bo y pi c 
is added in the 4-DOF case and the 12-DOF case includes the first eight elastic modes. 
Figure 4-61 shows the booster normal acceleration A and N 0 variations along the use- 
lage. Comparing these results with those for the booster return configuration shown in 
Figure 4-14 indicates that rms acceleration is generally higher during ascent but the 
characteristic frequency over the forward half of the booster is lower. Figures 4-62 
and 4-63 show the fuselage crew compartment vertical acceleration trans er nc ion 
and PSD plots. The A- matrix determinant magnitude and phase angle plots are shown 
in Figure 4-64 and in polar form in Figure 4-65. 

The load summary for the ascent flight, symmetric, SAS-on case is given in Table 4-7, 
•The rms loads are generaUy higher with SAS on than SAS off for the reasons previously 
mentioned. However, the loads are still well within the structural capabilities of the 
vehicle for design turbulence levels. For example, Figure 4-66 shows that to exceed 
the design limit values of wing root bending moment and shear, the rms turbulence 
velocity would have to be greater than 300 feet per second. The SAS-on crew station 
acceleration transfer function and PSD plots are contained in Figures 4-67 and 4-68, 
respectively. The A-matrix determinant is shown in Figures 4-69 and 4-70. 

Table 4-8 contains the load summary for the ascent flight, antisymmetric » SAS-off 
case. The 3-DOF case consists of the two slosh modes and rigid body lateral trans- 
lation. The 5-DOF case includes rigid body yaw and roll, and the 12-DOF condition 
includes the first seven elastic antisymmetric modes. The lateral acceleration A 
and N values for the booster fuselage are shown in Figure 4-71. Figures 4-72 and 
4-73 show the crew station lateral acceleration transfer function and PSD, respectively. 
The A-matrix determinant is shown in Figures 4-74 and 4-75. 

The SAS-on results for the antisymmetric ascent condition are summarized in Table 
4-9. The lateral acceleration transfer function and PSD for the crew station are 
shown in Figures 4-76 and 4-77, and the A-matrix determinant is shown in Figures 

4-78 and 4-79. 
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Table 4-6. Booster Ascent Flight Load Summary (SAS Off, Symmetric) 



3 DOF 


4 DOF 

12 DOF 

Response Item 

A 

N q (Hz) 

A 

N o (Hz) 

A 

N 0 (Hz) 

Wing Root Shear (lb) 

1,435 

0.800 

1,726 

0. 674 

1,680 

0.928 

Wing Root Bending Moment 
(in. -lb) 

0.407 X 10 6 

0. 846 

0.492 X 10 6 

0. 707 

0.477 X 10 6 

0.980 

Wing Root Torque About 

Station 3488. 5 (in. -lb) 

0. 110 X 10 6 

1. 545 

0. 164 X 10 6 

1. 033 

0.314 X 10 6 

3.864 

Fuselage Vertical Shear at 
Station 2800 (lb) 

1,225 

1. 923 

1,487 

1.612 

2,605 

2. 731 

Fuselage Vertical Bending 
Moment at Station 2800 (in. -lb) 

0. 76 X 10 6 

1. 956 

0. 729 X 10 6 

2.080 

2. 26 X 10 6 

2.660 

Forward Attachment Link 
Vertical Load (lb) 

18,515 

0. 866 

20,740 

0.767 

20,361 

0. 860 

Aft Attachment Link 
Vertical Load (lb) 

6,023 

1.283 

5,116 

1.519 

5,863 

2. 061 
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Figure 4-63. Ascent Flight Fuselage Normal Acceleration 
PSD. (SAS Off, Symmetric) 
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Figure 4-64. Ascent Flight A-Matrix Determinant Magnitude 
and Phase Angle (SAS Off, Symmetric) 
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Table 4-7. Booster Ascent Flight Load Summary (SAS On, Symmetric) 



3 DOF 

4 DOF 

12 DOF 

Response Item 

A 

N 0 (Hz) 

A 

N 0 (Hz) 

A 

N 0 (Hz) 

Wing Root Shear (lb) 

1,435 

0. 800 

1,486 

0. 783 

1,568 

1.316 

Wing Root Bending Moment 
(in. -lb) 

0.407 x 10 6 

0. 846 

0.417 X 10 6 

0. 831 

0.440 X 10 6 

1.350 

Wing Root Torque About 
Station 3488. 5 (in. -lb) 

0. 110 X 10 6 

1. 545 

0.548 

0. 704 

1. 623 X 10 6 

1. 783 ' ' 

Fuselage Vertical Shear 
at Station 2800 (lb) 

1,225 

1. 923 

7,955 

0. 708 

8,870 

1. 841 

Fuseiage Vertical Bending 
Moment at Station 2800 (in.-lb) 

0.760 X 10 6 

1.956 

4. 23 X 10 6 

0. 731 

3.924 X 10 6 

2. 348 

Forward Attachment Link 
Vertical Load (lb) 

18,515 

0. 866 

21,413 

0. 835 

32,193 

1.349 

Aft Attachment Link 
Vertical Load (lb) 

6,023 

1.283 

14,522 

0. 769 

15,039 

,1. 636 




FRACTION OF DESIGN LIMIT VALUE 



RMS TURBULENCE LEVEL (ft/sec) 

Figure 4-66. Ascent Flight Wing Root RMS Shear and Bending Moment 
Versus RMS Turbulence Level (SAS On, Symmetric) 
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Figure 4-67. Ascent Flight Fuselage Normal Acceleration 
Transfer Function (SAS On, Symmetric) 
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FOR 1.0 FT/SEC RMS GUST - G**2/CPS 



Figure 4-68. Ascent Flight Fuselage Normal Acceleration 
PSD (SAS On, Symmetric) 
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Table 4-8. Booster Ascent Flight Load Summary (SAS Off, Antisymmetric) 



3 DOF 

5 DOF 

12 DOF 

Response Item 

A 

N c (Hz) 

A 

N 0 (Hz) 

A 

N 0 (Hz) 

Wing Root Shear (lb) 

0 

0 

384 

0.313 

372.0 

1.858 

Wing Root Bending Moment 
(in. -lb) 

0 

0 

0. 137 X 10 6 

0.251 

0.114 X 10 6 

1. 598 

Wing Root Torque About 
Station 3488. 5 (in. -lb) 

0 

0 

32,308 

0. 702 

50,198 

2. 567 

Fuselage Lateral Shear 
at Station 2800 (lb) 

1,399 

1. 177 

3,096 

0.718 

3,378 

1. 166 

Fuselage Lateral Bending 
Moment at Station 2800 (in. -lb) 

1.67 X 10 6 

0.922 

‘ 

3.22 X 10 6 

0.812 

3. 59 X 10 6 

1.429 

Forward Attachment Link 
Lateral Load (lb) 

3,134 

1. 074 

2,908 

0. 849 

3,611 

1. 444 

Aft Attachment Link 
Lateral Load (lb) 

2,484 

1. 536 

653 

2.318 

1,470 . 

2. 118 




A (g/ft/sec) 



PHASE ANGLE ACCELERATION CGI PER 1.0 FT/SEC GUST, G/FPS 

COEGREES) 


PT. 2 ACC. AT F.S. 1150 CREW STATION 
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Figure 4-72. Ascent Flight Fuselage Lateral Acceleration 
Transfer Function (SAS Off, Antisymmetric) 
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Figure 4-73. Ascent Flight Fuselage Lateral Acceleration 
PSD (SAS Off, Antisymmetric) 
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Figure 4-74. Ascent Flight A-Matrix Determinant Magnitude 
and Phase Angle (SAS Off, Antisymmetric) 
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Table 4-9. Booster Ascent Flight Load Su mma ry (SAS On, Antisymmetric) 



3 DOF 

5 DOF 

12 DOF 

Response Item 

A 

N 0 (Hz) 

. A 

N q (Hz) 

A 

N 0 (Hz) 

Wing Root Shear (lb) 

0 

0 

2.40 

3. 712 

433 

2.985 

Wing Root Bending 
Moment (in. -lb) 

0 

0 

611 

3.651 

0. 109 X 10 6 

3. 000 

Wing Root Torque About 
Station 3488. 5 (in. -lb) 

0 

0 

• 

518 

3. 800 

57,600 

3. 980 

Fuselage Lateral Shear 
at Station 2800 (lb) 

1,399 

1. 177 

910 

0. 898 

1,315 

2.078 

Fuselage Lateral Bending 
Moment at Station 2800 (in. -lb) 

1.67 X 10 6 

0. 922 

1. 08 X 10 6 

0. 842 

1. 56 X 10 6 

2. 110 

Forward Attachment Link 
Lateral Load (lb) 

3,134 

1. 074 

2,475 

0. 899 

4,720 

3.341 

Aft Attachment Link 
Lateral Load (lb) 

2,484 

1.536 

3,233 

1. 167 

9,369 

3.358 
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Figure 4-76. Ascent Flight Fuselage Lateral Acceleration 
Transfer Function (SAS On, Antisymmetric) 








SECTION 5 


CONCLUSIONS AND RECOMMENDATIONS 


A computer program that can analyze the response of space shuttle to atmospheric 
turbulence has been developed and demonstrated. The method developed accounts for 
propellant slosh, gimballed engine, and stability augmentation system (SAS) coupling 
with the rigid body and elastic modes of the mated or separated vehicles. Statistical 
outputs relating vehicle loads and accelerations to the level of random turbulence are 
generated. Time histories due to discrete disturbances can also be obtained. 

By formulating the problem in the frequency domain, the capability is achieved to use 
any unsteady or quasi-steady aerodynamic theory based on harmonic motion for gen- 
erating aerodynamic forces due to both the turbulence and the vehicle response. Gust 
penetration effects are accounted for. The quasi-steady aerodynamics approach based 
on the method of Woodward appears to have good potential for the space shuttle turbu- 
lence response problem because of its ability to account for aerodynamic interference 
and body flow effects. It is recommended that Woodward theory be applied in future 
space shuttle turbulence response studies. 

Although load and acceleration responses calculated in the present analysis are well 
within vehicle design limits , the autopilot gains caused apparent instabilities that 
essentially invalidate the results of the SAS-on conditions. It is recommended that an 
elastic-vehicle/SAS-stability analysis preceed the turbulence response calculated in 
future studies to ensure that autopilot parameters are realistic. 

In the present analysis, only one ascent and one booster flyback flight condition were 
analyzed. This work should be extended to other flight times, concentrating on ascent 
flight where the complex aerodynamic and elastic problems are expected. 
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